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Abstract. - It is analytically shown that the one-dimensional Ising model with Glauber 
dynamics exhibits short time memory effects when submitted to an abrupt change in the 
temperature. These effects are qualitatively similar to those experimentally observed in the 
compaction of vibrated granular materials. Moreover, a critical time separating regimes of 
"normal" and "anomalous" responses to the perturbation is found. 



The dynamics of granular systems exhibits a rich phenomenology that has been the subject 
of intensive study in the last years jlj. In particular, when submitted to vertical vibration 
starting from a loose packing configuration, granular media compact behaving in a way that is 
reminiscent of conventional structural glasses The glassy nature of granular compaction 

clearly shows up when the response of the system to a sudden change in the vibration intensity 
is measured . The results show the presence of relevant short-term memory effects that are 
in contrast with the long time behaviour observed at constant intensity. A decrease (increase) 
in the vibration acceleration leads to an increase (decrease) of the compaction rate at short 
times. For large times, the previous history is forgotten and the relaxation rate becomes the 
same as for a constant intensity process. 

Here we study the existence of a similar "anomalous" response in the context of the 
one-dimensional Ising model with Glauber dynamics 0. This is a model with short range 
interactions that exhibits many of the qualitative features of conventional structural glasses, 
like slow non-exponential relaxation, a kinetic phenomenon resembling a laboratory glass 
transition, aging, and hysteresis effects |§-^9|. On the other hand, it is simple enough to 
allow exact analytical calculations in situations where more realistic models are untractable. 
In particular, we will be able to calculate the short-time response of the system to a sudden 
small change in the temperature, which plays in the model a role analogous to the vibration 
intensity in granular compaction. As long as the change is not made at the early stages of 
the relaxation process, the system shows an "anomalous" response. On the other hand, if 
the perturbation is introduced very soon, the response is "normal" , in the sense that the 
modification of the relaxation rate has the same sign as the temperature change. Let us 
mention that these two qualitatively different time regimes have been also identified within a 
general theoretical framework developed to analyze memory effects in granular compaction [|| , 
although it has not been experimentally verified yet to the best of our knowledge. 
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The energy of the one-dimensional Ising model in the configuration <x = {<7j}, Uj = ±1, is 

H((T) = -J^2(a k a k+1 -1), (1) 

k 

where J > is the ferromagnetic coupling constant between nearest neighbour spins. The 
Glauber dynamics for the model is formulated by means of a master equation, with the 
transition rate from configuration <x to configuration R k cr, differing from er only in the flip of 
spin a k , given by M 

^ T T 1 

Uk{&) = 1 _ ^k(crk-i +crfc+i) • (2) 

Here 

, 2J 

7 = tanh— — , (3) 
k b± 

ks being Boltzmann's constant and T the temperature of the system. The quantity a deter- 
mines the global time scale of the system. Although it can be taken as temperature indepen- 
dent, in the context of structural glasses, where spin configurations are associated with the 
minima of the potential energy Q4 |, an Arrhenius expression 

a{T) = a Q e^(—^\, (4) 

is often considered |l5|, . The constant B is a measure of the energy barrier separating 
energy minima in configuration space and ag can be employed to define a dimcnsionless basic 
time scale. In the following, we will use eq. (Q) for the sake of generality. The particular case 
of temperature independent a can be obtained by taking the limit B — > 0. 

We will restrict ourselves to homogeneous situations, i.e. states where all the averages are 
translationally invariant. It is useful to introduce the spin-spin correlations 

C n (t) = (a k a k+n ) t = ^ Vk<Jk+nP (a; t) , (5) 

where p(<r, t) is the probability of configuration er at time t. In an infinite lattice, the time 
evolution of these correlations at constant temperature is given by 

j t C n = -2aC n + ai {C n ^+C n+1 ), (6) 
for n > 1, while Cq = I. The equilibrium correlations are 

CW(T) = £(Tr, £( T )=tanh^. (7) 
Denoting by N the number of spins, a dimensionless energy per spin can be defined as 

E{t) = ^j = l-C 1 {t). (8) 

An equation for the time evolution of E(t) follows by particularizing eq. (Q) for n — 1. We 
will write it as an expression for the instantaneous energy relaxation rate r(t), 

r(t) = = a(T) M*) " e(T)v*(t)] , (9) 
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where 

fi 1 (t)=2E(t)-l + C 2 (t)>0, /z 2 (f) = 1 + C 2 {t) > 0, (10) 
and we have introduced 

e(T) = l- 7 (T). (11) 

Both e(T) and a(T) are monotonic increasing functions of the temperature. An equation 
having the same structure as eq. (|^) was proposed to describe compaction in tapped granular 
systems in ref. The main difference is that here we have considered the time evolution 
of the energy, instead of the density. In the interpretation of the Ising model as describing 
structural relaxation, an increase of the density corresponds to a decrease of the energy. This 
is the reason why we have defined r(t) with a minus sign in eq. (^|). 

Let us suppose we carry out the following experiment. Starting from a given configuration, 
the system is allowed to relax at constant temperature T. At a certain time t w , before the 
system has reached the steady configuration corresponding to T, the temperature is suddenly 
changed to T + AT. This produces a jump Ar w in the relaxation rate. For AT <C T, keeping 
only first order terms in the perturbation AT, it is easily obtained that 

Ar w = X(t w )Aa, (12) 

with 

r(t) , \ de r(t) nil - ■y 2 ) , . 

a da a I 

and r\ = AJ/B. The function X(t) is defined along the relaxation curve at constant T. As long 
as AT is small enough so that the linear approximation we are considering is accurate, the 
function X(t) for t = t w determines the relative behaviour of the relaxation rate jump Ar w 
with respect to Aa (or AT). 

To study the function X(t) we have to evaluate the relaxation rate r(t) and the correlation 
function /^(i) for a process at constant temperature. The general solution of the hierarchy (||) 
for an arbitrary initial condition can be analytically found [0,^l|. We will consider processes 
starting in the state of maximum disorder, i.e. the equilibrium state at infinite temperature. 
Similar results are obtained for other initial conditions, as long as they correspond to highly 
disordered states. A simple calculation gives [Q,|2l| 

A n (t)^C n (t)-C^(T) = -l r rfg ^gsinn g e _ 2at(1 _ 7COBg) 

TT JO l-7 cos <7 

Use of the above expression with n = 1 into eq. M) yields 

E(t)-E^(T) = 1 dq S ' n2<? e -2°*(i-7°°»g) (15) 
7r / 1 - 7 cos q 

where (T) = 1 — c[ (T) is the equilibrium value of the dimensionless energy per particle. 
The relaxation rate of the energy, defined in eq. ([I), for the present process is 



r(t) — AAL I dq sm 2 qe- 2at(1 -< cos ^ = A^AAAl > ( i6 ) 



2c f7 r^^2„„-2a t (l-^cos«) e- 2at h{2 al t) 

n 

where I\ is the modified Bessel function of the first kind . 

To evaluate n^it) we use again eq. (|l4]), now for n = 2. In this way we get from eq. ( |l0| ) 

^(t) = i+e- 2 f(t), (17) 

7T 
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dg sing sin 2g e _ 2at(1 _ 7COBg) (lg) 
o 1 - 7 cos q 



where 



Substitution of eqs. ( J16| ) and ( |17| ) into eq. (|13|) leads to 

at 2 L 7r J 

In the following we will consider that the temperature T is very low, so that s <C 1, 7 — ► 1, 
£ — ► 1, and a — > 0. It could be concluded that the second term on the right hand side of 
eq. ([[9]) can be neglected as compared with the first one |22| in this limit . Nevertheless, this 
is not true in general, and the relative relevance of both terms depends on the time scale of 
interest. In particular, for t — > 00 the first term vanishes, as well as f(t), and 

Xoo = lim A(t) = -- 1 ' { — ; ~ -2 V s < 0. (20) 

t — >OQ 2 

Since X(t) is a continuous function of t, the above result implies that, for large enough t Wl 
X{t w ) < and, therefore, Ar w and Aq have opposite signs. A sudden increase of the tem- 
perature produces a negative jump of the relaxation rate. This is the kind of "anomalous" 
response found in compaction processes of granular materials . 
It is convenient to define the function 

at Ztt 

that vanishes for t — ► 00. We are going to analyze the behaviour of this quantity both in the 
short and the large time regions. First, consider "short" times characterized by 

at- O(l), (22) 

so that sat <C 1 for the low temperatures we are dealing with. On the time scale defined 
above, we can approximate 7 by unity in the expression of f(t), eq. (p"8|), obtaining 

f(t) = ne~ 2at [h{2at) + 2h{2at) + 7 2 (2crt)] . (23) 

When this expression is substituted into eq. ( pl| ) , it is easily realized that the term involving 
f{t) can be neglected because of the r\e factor. Then, 

^ f^iiM, (24) 

at 

for at — O(l), £ < 1. In this time window, X(t) — A(t) > 0, where we have taken into 



account that Aoo is of order e (see eq. (20)). As a consequence, the response of the system to 
a small change AT in the temperature at t = t w is "normal" if aet w <C 1. The modification 
of the relaxation rate has the same sign as the temperature change. Above, it was shown that 
for asymptotically large times the response was qualitatively different. Let us study in more 
detail what happens in the "large" time region. We define a slow time scale by 

T = aet~ 0(1), (25) 

i.e., times of the order of the system relaxation time. For these times, it is at — Q(e^ 1 ) 3> 1 
and we can approximate Ii(2ajt) by its asymptotic behaviour for large argument |23fl , getting 

e- 2at h(2ont) ,3/2 e~ 2 - 

at 27rV2 T 3/2- W 
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This quantity is much smaller than Aoo. Moreover, a simple asymptotic analysis gives 

/(0^(2^) 1 / 2 r(-l/2,2r), (27) 

r(a, x) being the incomplete Gamma function p3]. Use of eqs. (|26| ) and ( p7| ) into cq. 
leads to 



A(i) 



.3/2 



-2t 



27r l/2 T 3/2 



1/2 



r(-l/2,2r) 



(28) 



valid for e <C 1, r = O(l). In this long time regime, A(i) is of the order of e 3 / 2 , and 
|Aoo| » A(i). Then, X(t) ~ < 0. The response of the system is anomalous, as we already 
knew from the discussion of the asymptotic behaviour for t — > oo. 

Since we have found that X(t) has different signs in the short and long time regions, because 
of continuity reasons there must be at least one critical time t c at which X(t c ) = 0. This time 
separates regions with normal and anomalous behaviour. Furthermore, the above discussion 
indicates that this time must belong to an intermediate time scale between those previously 
considered, namely those given by eqs. (p2| ) and (p5|), respectively. The matching time range 
where the two solutions, eqs. (E4|) and (|28|), are valid is characterized by 



at > 1, r = eat <C 1. 

On this time scale, we can substitute the Bcssel function in eq. (^ 
for large argument J23|, with the result 



A(t) 



1 



27T 1 /2( Q , f )3/2 



(29) 

by its asymptotic behaviour 



(30) 



The same expression is obtained from eq. (|2q ) in the appropriate limit, showing that the 
solutions corresponding to the short and large time scales match in the intermediate scale 
dchned by eq. ( |29|) . A single uniform approximation A un if(t), valid for all times, is obtained 
by adding eqs. ^2J) and (p8|), and subtracting the expression in the matching region, eq. (|30|), 



A-unif (^) 



t h{2a£) 



^ — 2ast 



- 1 



at 



2w 1 / 2 (at) 3 / 2 



e 3 ^ i - 



1/2 



r(-l/2,2aet). 



(31) 



In fig. [l] we compare this result with the exact expression, eq. ( |2l| ) for s = 0.01 and J = B, 
i.e. r\ = 4. A very good agreement is observed for all times. The relevant physical quantity 
we are interested in is the critical time t c , verifying X(t c ) — or, equivalently, A(i c ) = —Aoo. 
This latter quantity is also indicated in the figure. Since t c is expected to belong to the 
intermediate scale, we get from eqs. (p0|) and (30) 



at c 



-2/3 



(32) 



This result confirms that t c is in the overlapping region in agreement with our previous ansatz, 
since in the limit e — > it is at c 3> 1 and r c = eat c <C 1. Therefore, eq. (133) is consistent 



with the uniform expression (31). It is interesting to analyze the temperature dependence of 
t c . By taking logarithms in cq. (E32j) we find for low T, 



ln(a *c 



■In 



B 
~k~RT 



2rj 

y 



(33) 
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Fig. 1 Fig. 2 

Fig. 1 - Comparison of the numerical evaluation of eq. (j2l| (circles) with the uniform asymptotic 
approximation given by eq. ( fjl] ) (solid line). Also plotted (dotted line) is the asymptotic long time 
value |Aoo|- Note that the quantities represented in both axes are dimensionless. 

Fig. 2 - The dimensionless critical time t c , separating the regions of normal and anomalous response, 
as a function of the inverse of the reduced temperature. The circles have been obtained by numerical 
evaluation of eq. (pl|), while the solid line is the result of the asymptotic calculation, eq. (^). 



showing an Arrhenius dependence. In fig. || the logarithm of the critical time t c , given by 
eq. (|||), is plotted as a function of the inverse of the reduced temperature B/ksT, again for 
i] = 4. For comparison, also the numerical values obtained from the exact expression, eq. (^l|), 
are included. A linear profile, consistent with eq. (|33|), is observed. 

Let us now address the particular case of a being temperature independent that can be 
obtained by considering the limit B — > 0. Then r\ goes to infinity and t c vanishes. Therefore, 
the system always exhibits an anomalous response to a temperature jump, independently of 
the time instant in which it is produced. The physical reason for this result is clear, since eq. 
( fl2| ) reduces in this situation to 

Ar w de 

-£=■ = -a/i 2li ,— < 0. (34) 

Equation ( |3~i| ) reflects that, if a does not depend on the temperature, the rates of the processes 
decreasing the energy do not change instantaneously in a temperature jump. This can be 
directly verified, for instance, from eq. (0). 

In conclusion, the major result obtained in this paper is the analytical derivation of memory 
effects in the relaxation of the one-dimensional Ising model with Glauber dynamics. More- 
over, these effects are qualitatively similar to those observed experimentally in granular com- 
paction 0-^, for which the existence of a critical time separating the regions of "normal" 
and "anomalous" responses, has been theoretically argued Q. Therefore, also in the present 
context the Ising model seems to be the simplest relevant system beyond mean-field descrip- 
tions. The results in this paper seems to indicate that the memory effects and the presence of 
a critical time, as discussed here, may be quite general. In the Ising model, the critical time 
lies in an intermediate time window, being much smaller than the average relaxation time of 
the system towards its equilibrium state. 
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